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ABSTRACT 

a 

We report calculation of critical exponents corresponding to three thermodynamic func- 

CT, 



tions, viz. chemical potential, condensate fraction and specific heat of a non-interacting 
Bose gas confined in an isotropic harmonic oscillator trap in one-, two- and three- 
dimensions. Critical exponents are also calculated for bosons interacting through van 
der Waals potential and confined in a three-dimensional harmonic trap. These expo- 
nents are found to be independent of the size of the system and whether the bosons 

£> ' are interacting or not, demonstrating their universality property. We also find that a 

VO - 

CN . continuous phase transition occurs even in one-dimension, as shown by a rapid change 

0\ ■ 

of the thermodynamic functions at the critical point. This is in sharp contrast with 

q ■ the result found in the text book treatment of replacing the sum over occupations of 



O 



X 



discrete states by an integral. The two- and three-dimensional results agree with the well 
known results in the thermodynamic limit, i.e. when the number of bosons in the trap 
approaches infinity. In the same limit, the criticality disappears for the one-dimensional 
trap, which again agrees with the standard result. 

PACS numbers : 03.75.Hh, 03.75.Nt, 64.60.F-, 75.40. Cx 

Bose-Einstein condensation (BEC) is the transition process, in which a macro- 
scopic fraction of bosons goes into the lowest energy state, below a certain critical tem- 
perature (T c ) [l|. It was predicted by Einstein in 1925, based on Bose's explanation 



of black body radiation. The BEC is generally considered as a phase transition in the 
thermodynamic limit. In standard text books |l|, |2j, it is shown that BEC is possible 
in an ideal uniform (i.e. not confined) Bose gas at a finite temperature in three and 
higher dimensions, while in two-dimension the transition temperature becomes zero and 
no BEC transition is possible in one- dimension. In a uniform ideal (non-interacting) Bose 
gas in three-dimension, the condensate fraction (defined as j$-, where N is the number of 
bosons in the lowest energy state and N is the total number of bosons) at a temperature 
T < T c is given by jj 



— = l-(-) § (1) 



while the specific heat is given by 



m 



C N {T) = -Nk B )^-(-)->, (2) 



4 aiyT c 

where ks is the Boltzmann constant and ((x) = Y^=i n ~ x * s the Riemann zeta function. 
Both condensate fraction and Cn are continuous, while the first derivatives are discon- 
tinuous at T = T c for the ideal Bose gas |l|. This means that in the thermodynamic 
limit (N —?• oo) BEC is a continuous phase transition process. The situation is quite 
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larmonic oscillator traps. In 
Thus, the corresponding 



different if the Bose gas is non- uniform, i.e., confined in a suitable trap [4]. Commonly 
used experimental traps are spherical or axially symmetric 
such cases, Cn(T) has a finite discontinuity at T = T c |l|, |4j. 
phase transition is of first order. However, it is important to note that in the usual 
text book treatment, the sum representing occupation of single particle levels is replaced 
by an integral, involving the density of states which is obtained by a semi-classical ap- 
proximation. In doing so, the occupation of the lowest energy level is taken separately. 
This approximation is valid only in the thermodynamic limit of an ideal uniform Bose 
gas (having continuous energy levels). In a potential trap, with discrete energy levels, 
this replacement is not justified for a finite number of atoms. Replacing the sum by an 
integration including semi-classical density of states, Bagnato and Kleppner [5( showed 
that BEC is possible in one- dimension, if the confining potential is more confining than a 
parabolic potential, while in two- dimension, BEC is possible for anypower law confining 
potential. Using a modified semi-classical approximation, Yukalov |6| replaced the lower 
limit of integration from zero to a finite value given by the uncertainty relation. He thus 



removed the divergence at the lower limit of the integration to obtain Bose condensation 
for an arbitrary power law confinement in any dimension. For a finite value of N, in con- 
trast with the text book treatment of an ideal Bose gas, Ketterle and van Druten treated 
the sum directly and showed that BEC is possible even in one-dimension, the critical 
temperature being higher and therefore more convenient for experimental achievement 
of BEC in quasi-one dimensional traps Ml. In such cases, both the condensate fraction 
and the specific heat are smooth and continuous functions of T near T c [7|. In two and 
three dimensions, the specific heat at constant particle number and for a fixed confining 
potential attains its maximum at the transition temperature (T c ) and decreases rapidly 
thereafter. For the three dimensional case, its value approaches the classical value fiVfee 
(for ideal uniform Bose gas) or 3Nks (for non-interacting atoms trapped in a 3D har- 
monic oscillator well). For large N, their behavior approaches that in the thermodynamic 



limit 



a. 



Other studies also showed that BEC in lower dimensions can be obtained using 



a thermodynamic limit approach 



B,Q. 



An important physical quantity in connection with the phase transition at a 

critical point is the critical exponent. A critical point is characterized by divergences, 

discontinuities or rapid changes in thermodynamic quantities, depending on the nature 

of the thermodynamic quantity, dimensionality and the size of the system. The critical 

exponent is useful in understanding the rapidly changing behavior of thermodynamic 

functions at the critical temperature. In terms of a reduced temperature 

T — T 
«=^, (3) 

which is a dimensionless measure of the deviation of temperature from the critical tem- 
perature, one can define the critical exponent (A) for a thermodynamic function F(t) 



as [11 1 



A = lim™. (4) 

t->o ln|t| 

This corresponds to F(t) ~ |t| A , for small values of \t\. Obviously, at the critical temper- 
ature, F(t) must either vanish or be singular, depending on the sign of A. For quantities 
that are neither singular, nor vanish at T c , (e.g., specific heat, condensate fraction or 
chemical potential of a BEC containing a finite number of atoms), one can take F(t) to 



be the difference of the desired thermodynamic function at T and its value at T c . Indeed, 
in this case, A is zero and our prescription will give Ai of the relation 

F(t) = F(0) + bt Xl + ... , (5) 



for positive Ai (see Eq. (2.21) of Ref. [ill]). Since in this case, A is zero and Ai is the 
leading exponent of interest, the nature of the thermodynamic function near the critical 
point is determined by Ai. Hence for such a situation, we define X 1 as the critical exponent 
(A). Therefore our modified definition is 

A=Bni H£ev£pai. (6) 

t->o ln|t| 

An important characteristic property is the universality of the critical exponent. 
The value of A does not depend on the interatomic interaction or detailed nature of the 
system. Its value dep ends only on the dimension of the system and the symmetry of the 



order parameter llj. In the present work, we calculate critical exponents for a number 
of thermodynamic functions of the BEC in different cases and investigate whether they 
depend on the system parameters. We consider a non-interacting Bose gas trapped 
in a harmonic oscillator potential in one, two and three dimensions. Finally, we also 
include realistic interatomic interactions in an approximate many-body treatment of 
the real Bose gas. As discussed earlier, the thermodynamic functions like specific heat, 
condensate fraction and chemical potential depend on the dimension of the system, the 
choice of the trap potential, whether interatomic interactions are included or not, etc. 
For a proper phase transition, it is expected that the critical exponents will show their 
universal characteristics. Our present work explores this possibility. 

We consider N non-interacting bosonic atoms trapped in an isotropic harmonic 
potential of frequency uj. The energy scale is so chosen that the single particle ground 
state is at zero energy. The energy eigenvalues E n (n = 0, 1, 2, ...) are given by 

E n = nhu (7) 

The average number of particles in each single particle energy eigenstate at an absolute 



temperature T i.e., the Bose distribution function f(E n ) is given by 

^») = e /?(^-M) _ 1 ^ 

where (3 = l/ksT and \x is the chemical potential. The latter is determined from the 
constraint that the total number of particles is N. Clearly, jjl has a temperature de- 
pendence. In general for a d- dimensional space, an energy level is degenerate with a 
degeneracy factor j n . The degeneracy factors for the n-th level of the one-, two- and 
three-dimensional harmonic oscillator are 1, (n+1) and — — -^ — - respectively. The total 
number of bosons in the trap is fixed and can be written as 

+oo 

N = J2lnf(E n ) (9) 

n=0 

The total energy for the system is given by 

E(N,T) = Y,lnf(E n )E n (10) 

ra=0 

The specific heat at fixed particle number (N) is calculated using the relation 

c H m = B2 (ii) 

Using (jSJ), ffTUl) . (fTTj) one can obtain the heat capacity as 

r* t r r\ — a v^+°° inE n exp {p{E n -ii)) \ E n -\i ,9^1 n o\ 



n=0 (exp(^(£„-/i))-l) 2 L T ' 8T 

+oo 7n E n exp (p(E n -ii)) , E n -n _ E^ 7m(gm-M)exp(/3(E m - At ))(/(E m )) 2 
n=0 (exp(^(B„-^))-l)2L T TE+^ 7p exp ( / 3(E J) - /U ))(/(B p ))2 



For an ideal bosonic gas containing iV bosons in a three-dimensional isotropic 
harmonic well, there is a well defined critical temperature T c ° (which can be treated as a 
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reference temperature in what follows) and is given by [3] 

where £(3) ~ 1.202057. This is obtained from the condition that all the particles are 
in the lowest single particle state for T < T 6 °. Hence \x remains equal to the energy of 
the single particle ground state for T < T c ° and start decreasing rapidly for T > T c °. 
However as mentioned earlier, for a correct treatment [i.e., not replacing the sum in 
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by an integral], the number of bosons in the ground state becomes non-zero for 
any temperature above zero. Consequently, \x decreases slowly from its maximum value 
(equal to the ground state energy) as T increases from zero, the rate of increase becoming 
suddenly rapid at some temperature close to the reference temperature T c ° corresponding 
to the same value of N. Thus, in this case there is no sharp critical temperature. The heat 
capacity also becomes a smooth function of T, attaining a maximum at a temperature, 
at which \x suddenly becomes a rapidly decreasing function of T (except for the one- 
dimensional case, see below). The critical temperature T c is defined as the temperature 
at which Cn(T) is a maximum [8( 

9Cn(T) _ 

For a finite value of N, T c is less than T c ° and approaches T c ° as N — > oo. This is 

expected, since in this limit, the approximation of replacing the sum in (Q by an integral 

becomes exact. Using our adopted definition, the critical exponent for the heat capacity 

is given by 

A s P ht = Um \n\C N (t)-CT x )\ (15) 

t-+o- In |t| 

where C™ ax is the maximum of Cjv (T) and t is given by Eq. (J3]) • 

In a similar fashion, the critical exponents for the chemical potential and the 
condensate fraction is obtained. For all thermodynamic functions, we use the same 
critical temperature defined by (JHJ). Since the thermodynamic functions considered here 
are all continuous at T = T c , we take the difference of the value of the function at T with 
that at T c . The critical exponent for the chemical potential is thus given by 



A chem P ot = Um In |/i(T) - //(T t 
and that for the condensate fraction is 



, . . (16) 

t^o+ In t K ' 



A Condfrac = ^ \ N N 1 _ (1?) 

i^o- In |t| 

For the interacting case, one has to solve the iV-body Schrodinger equation. 
An exact solution is an almost impossible task. The mean-field theory together with 
the assumption of a contact two-body interaction leads to the commonly used Gross- 



Pitaevskii equation (GPE) 



12| . In this approach, all correlations are disregarded and 



no realistic two-body interaction can be used. In addition, for an attractive BEC, there 
is a pathological s in g n.a„t y at the origin Q Instead, we adopt amplified few-body 
technique, called potential harmonics expansion method (PHEM) |14l. |15|. in which all 
two-body correlations are retained (disregard of higher-order correlations is manifestly 
justified in a dilute BEC) and any realistic interatomic interaction can be incorporated. 
This technique has been successfully applied to both repulsive and attractive BECs [161 ]. 
In the following, we review the PHEM very briefly. Interested readers can get 
the details from Refs. J14J . Il5| . In this technique, center of mass motion is separated by 
introducing Jacobi coordinates, and the resulting Schrodinger equation for the relative 
motion is rewritten in terms of hyperspherical variables: the hyperradius (r) and (3iV — 4) 
hyperangles, collectively denoted by fijv-i [13] ■ Expansion of the relative wave function 



in the complete set of hyperspherical harmonics (HH) 17| leads to an exact treatment 



but involves (3N — 3) variables, which remains an almost impossible task for iV > 3. We 



choose the subset of HH, called potential harmonics (PH) subset [18| , which is sufficient 
for the expansion of the two-body interaction V(rij) in the HH basis, where f^ is the 
separation vector from the i-th to the j-th particle. Since the probability of three or more 
particles to come close to each other is practically zero in a typical BEC (indeed this 
is assured in an experimental BEC, in order to avoid three-body collisions, which lead 
to molecule formation and consequent depletion of the BEC), correlations higher than 
two-body ones are practically absent. Thus the expansion of each Faddeev component 
of the wave function in the PH basis, which retains two-body correlations, is manifestly 
appropriate. Substitution of this expansion in the relative Schrodinger equation and 
projection on a particular PH leads to a system of coupled differential equations in r. 
The latter is solved by hyperspherical adiabatic approximation [l9|, to obtain a large 
number of energy eigenvalues of the entire BEC for different orbital angular momenta 
of the system. These are then used in Eq.(j9]) to calculate the chemical potential at a 
chosen temperature T. Note that in this case, the sums over n in Eqs,©, f JTOj) and (12) 
are to be replaced by double sums over {n, 1} and 7„ is replaced by (2/ + 1). The energy 
eigenvalue E n \ is the energy of the n-th excitation for the orbital angular momentum / 



of the system. Finally Eq.(12) is used to calculate Cn{T) [20]. Condensate fraction is 



obtained as ^, where N is the n = term of Eq.(j9]). 



Initially, we consider non-interacting bosons in an isotropic harmonic trap in 
one, two and three dimensions. We first present the results for N non-interacting trapped 
bosons in a one- dimensional harmonic oscillator potential. In panel (a) of Fig. 1, we plot 
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FIG. 1. (Color online) Plot of chemical potential (panel a), condensate fraction (panel b) and 
Cn /Nks (panel c) as functions of ksl '/Hcu for indicated number (N) of non-interacting bosons 
in a one-dimensional harmonic trap. Note that chemical potential and ksT are expressed in 
energy oscillator unit (Hlj). Note also that Cn/NUb is a monotonically increasing function of 
T, giving rise to the common notion that there is no criticality in one dimension. 

fi/tiu against k B T/hio for iV = 1000, 10000 and 100000. For N = 1000, one notices that 
H remains nearly equal to its value at T = up to ksT/hw ~ 200, after which /i starts 
decreasing rapidly. This indicates a critical behavior. For larger number of particles, 
such a rapid change occurs at a much higher temperature, indicating that the critical 
temperature increases rapidly with N. This observation agrees with Ref. |3|. In panel 
(b) of Fig. 1, we plot the condensate fraction, as a function of ksT/hu for the chosen 
values of N. Once again, we notice that jS- has a rapid change at ksT '/ftw ~ 200 for 
iV = 1000, while for larger N the rapid change occurs at a much higher temperature. All 
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these indicate that there is a continuous phase transition at the critical point. However, 
from a plot of Cjy/Nks against k B T '/fvjj in panel (c) of Fig. 1, we notice that for the 
one- dimensional case, C^(T) is a monotonic function of T and so definition ( |14p cannot 
be used for calculation of T c . This implies that the one-dimensional case is distinctly 
different from higher dimensions (described later). This is reflected in the text book 
treatment (where sums are replaced by integrals), for which T c ° diverges since £(1) di- 
verges (compare with Eq. ( lT3"j) and Ref. J3(). Instead, we also notice from Fig. 1(c), that 
Cn{T) is almost linear for T far away from the critical region with distinctly different 
slopes. Hence for the one-dimensional case only, we define T c as the temperature at 
which the extrapolated asymptotic straight lines intercept. This gives T c = 243.77 nK 
for N = 1000. For N = 10000 and 100000 T c is larger (2132.5 nK and 15888 nK 100000 
respectively). The slopes of the Ctq{T)/NkB versus k B T /ftw curve before and after the 
critical point become closer to each other, as the particle number increases. Conse- 
quently, criticality gets smeared out as iV increases. Although existence of a continuous 
phase transition for finite N in a one-dimensional harmonic oscillator trap is in sharp 
contrast with common acceptance, the fact that the criticality gradually fades away as 
iV — > oo is in agreement. Using this value of T c for a particular value of N, we plot 
In \fi(T) - fi{T c ) |, In | ^^-^^ | and In \(C N (T)-C N (T c ))/Nk B \ against In \t\ in panels 
(a), (b) and (c) of Fig. 2 respectively, for the one dimensional trap. As expected, we 
obtain straight line fits. All such straight line fits for different N are parallel (to within 
calculation errors), showing that the critical exponents are independent of N. A typical 
straight line fit for each plot is shown by a black dotted line. The slopes of these lines 
for panels (a), (b) and (c) are respectively 0.919,1.122 and 1.004. These are then the 
critical exponents for chemical potential, condensate fraction and specific heat in the one 
dimensional harmonic oscillator trap. 

We next repeat the calculations for a two-dimensional isotropic harmonic 
oscillator trap. In this case, Cn{T) versus T curve has a maximum, from which we 
calculate T c according to Eq. (TT4"1) . Panels (a), (b), (c) of Fig. 3 present plots of 
ln|(/i(T) - fi(T c ))/hw\, \n\N (T)/N - N (T C )/N\ and ln\(C N (T) - C N (T c ))/Nk B \ re- 
spectively as functions of In \t\ for N = 1000, 10000, 100000 and 1000000 atoms in the 



10 







2 


N=10000 \~,*«- 
N= 100000 ---s^**? 
0.91 9*x+2. 092 - - / 


1 


/ ■■' 




//j^ 0m 





/ M 




/ * 


1 


" / / 




/ * 

/ i .• i 







-2 

ln|t| 

(a) 











N=10000 
N=100000 * J 
1.122*x-1.100 # 


wJ-/ 






- 




- 


-4 -2 






ln|t| 
(b) 




FIG. 2. (Color online) Plot of In |(/i(T)-/x(T c ))/M (P anel a )> ln \N (T)/N-N (T C )/N\ (panel 
b) and ln|(Civ(r) — CN{T C ))/Nks\ (panel c) as functions of ln|i| for indicated number (N) 
of non-interacting bosons in a one-dimensional harmonic trap. Typical straight line fits to 
the asymptotic regions are shown by black dotted lines. Parameter values of the fits are also 
indicated. Note that the asymptotic linear regions for different iV are parallel, showing that 
the critical exponent of a particular thermodynamic function is independent of N. 

trap. The asymptotic part of the plots are straight lines, which are parallel for different 
N of the same panel. This again indicates that the critical exponents are independent 
of N (note that a few points for large negative ln \t\ in panel (c) do not fall on the 
straight lines due to numerical errors, and have been disregarded). Typical straight line 
fits are indicated by black dotted lines and the fitting parameters are also displayed. 
From these, we find that the critical exponents for chemical potential, condensate frac- 
tion and specific heat are 0.915, 1.006 and 1.810 respectively for the two dimensional trap. 



We next calculate these quantities for a three dimensional harmonic oscillator 
trap. In Fig. 4, we plot fj,/Hu>, condensate fraction and Cn /Nks as functions of T/T® 
in panels (a), (b) and (c) respectively, where T c ° is calculated from Eq. f|T3|) . Since in 
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FIG. 3. (Color online) Plot of In |0i(T)-/x(T c ))/M (P anel a ), ln \N (T)/N-N (T C )/N\ (panel 
b) and ln \(Cn(T) — Ctf(T c ))/NkB\ (panel c) as functions of ln |i| for indicated number (JV) of 
non-interacting bosons in a two-dimensional harmonic trap. For other comments, see caption 
of Fig. 2. 



this case, Cn(T) shows a clear maximum, we can easily calculate corresponding T c from 
Eq. flU). Calculated values of TJT® are 0.896, 0.949, 0.976, and 0.989 for iV = 1000, 
10000, 100000 and 1000000 respectively. Using these values of T c , we calculate the differ- 
ences of the three thermodynamic functions from their values at the critical points and 
plot the logarithm of their moduli against ln \t\ in panels (a), (b) and (c) of Fig. 5. All 
calculated points for the chemical potential fall nicely on straight lines in the interval 
— 10 < ln |t| < —6. Interestingly, all curves for different N fall on the same straight line in 
this ln \t\ interval. This shows that the corresponding critical exponent is independent of 
N. For ln \t\ increasing above —6, the curves for different N separate gradually, showing 
that the nature of \x versus T curves away from the critical point depend on N. On the 
other hand, numerical errors start to show up for ln \t\ < —10. Hence such points are 
ignored for calculation of the critical exponent. The same behavior is seen for condensate 
fraction [Fig. 5(b)] and specific heat [Fig. 5(c)]. However the numerical errors are larger 
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FIG. 4. (Color online) Plot of chemical potential (panel a), condensate fraction (panel b) and 
Cn /NUb (panel c) as functions of T/T® [where T® is defined in Eq. (|13p ] for indicated number 
(N) of non-interacting bosons in a three-dimensional isotropic harmonic trap. Note that the 
chemical potential is expressed in energy oscillator unit (hui). In contrast with Fig. 1(c), panel 
(c) above shows a maximum followed by a sharp drop in Cjv, the sharpness of the drop increasing 
with N. Furthermore, as N — > oo, the temperature at which this sharp drop occurs approaches 

A c ■ 

for the calculation of the specific heat, showing some divergences. In the last case, the 
linear parts for different iV are nearly parallel showing that the corresponding critical 
exponent is independent of N. Typical straight line fits are shown by black dotted 
lines in Fig. 5, together with the fitted parameter values. Calculated critical exponents 
for chemical potential, condensate fraction and specific heat are 0.913, 0.995 and 1.900 
respectively. These values are independent of N, confirming universality of the critical 
exponent. 



Next we present the results of our calculation for iV bosons interacting via the 
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FIG. 5. (Color online) Plot of In \(p(T)-n(T e ))/hu\ (panel a), In \N (T)/N-N (T C )/N\ (panel 
b) and In \(Cn(T) — C]y(T c ))/NkB\ (panel c) as functions of In |i| for indicated number (N) of 
non-interacting bosons in a three-dimensional harmonic trap. For other comments, see caption 
of Fig. 2. 



van der Waals potential and trapped in an isotropic three dimensional harmonic oscilla- 
tor potential. As mentioned earlier, the GPE uses a contact interaction, whose strength 
depends only on the s-wave scattering length a s and hence is shape independent of the 
two-body potential. An earlier calculation [2l| showed that calculated observables are 
indeed shape dependent. Hence it is appropriate that a finite-range realistic interatomic 
potential like the vdW potential should be used. Calculated values of In |/i(T) — /i(T c )|, 

In | °^ ' iv I an< ^ ^ n l ( -^ Af (-^) — Cn(T c )\ J20(] have been plotted against In \t\ in panels 

(a), (b) and (c) of Fig. 6. Once again these plots are straight lines in the asymptotic 
(In \t\ — > — oo) region. However, due to large error in the numerical calculation, especially 
when \t\ is very small there are some fluctuations for large negative ln|t|. Hence the 
lower limit of In \t\ has been restricted to —8. From the linear plots we get the values 
of A Chempot = 0.913, A Condfrac = 1.005 and A Spht = 1.878. We notice that these values 
are independent of iV (since the asymptotic linear parts are parallel for different N, 
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FIG. 6. (Color online) Plot of ln|(/x(T) - fi(T e ))/fiu\ (panel a), ln\N (T)/N - N (T C )/N\ 
(panel b) and In |(Cjv(T) — CN(T c ))/NkB\ (panel c) as functions of In \t\ for indicated number 
of atoms interacting via van der Waals potential in a three-dimensional harmonic trap. For 
other comments, see caption of Fig. 2. 

in each of the panels (a), (b) and (c) of Fig. 6. We also notice that the values of the 
critical exponents for the interacting three dimensional case are equal (within numerical 
errors) to the corresponding critical exponents for the non-interacting three dimensional 
case. These observations show that the critical exponents of the BEC obey universality 
property. 



In conclusion, we remark that a continuous phase transition occurs in a non- 
interacting Bose gas trapped in a harmonic oscillator potential even in one- dimension. 
This is in sharp contrast with standard text book results [l|, |3|], but in agreement with 
Ketterle and van Druten |7|]. This is because the replacement of sums over single particle 
states by an integral over energy is not a valid approximation for discrete energy levels, 
especially when N is small. However, we find that there are some distinct characteristic 
features in the one-dimensional case. For example, the specific heat is a monotonic 
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increasing function, whereas in higher dimensions it has a maximum. The dependence 
of chemical potential, condensate fraction and specific heat on temperature for the non- 
interacting one-dimensionally trapped Bose gas has similar behavior as when a Bose 
gas interacts through the harmonic Calogero interaction [22|]. We prescribe a way to 
calculate the transition temperature (T c ) in this case (as the temperature at which the 
asymptotically linear curves for T far below and far above the critical region intercept). 
As N increases, T c increases rapidly. For N — > oo, the two asymptotic straight lines 
have the same slope and no T c can be defined. This agrees with the commonly accepted 
idea that there is no Bose-Einstein condensation in one- dimension. Ketterle and van 
Druten argued that a quasi one-dimensional experimental setup is advantageous, since 
T c becomes larger |Zj. However, in view of the above discussion, it is clear that such an 
experimental setup is useful for relatively small N. Although T c increases rapidly with 
N, the difference in the slopes of the two asymptotic regions of Cn{T) versus T curve 
decreases. Hence the difference in the properties of the two phases gradually reduce, 
making a distinct phase transition obscure. 

We have also calculated the critical exponents for three thermodynamic functions, 
viz. chemical potential, condensate fraction and specific heat at constant particle num- 
ber for a non-interacting Bose gas in an isotropic harmonic trap in one, two and three 
dimensions. For the three dimensional case, we have also included a realistic two-body 
interaction. We find that the critical exponent for a particular thermodynamic function 
does not depend on iV or whether the bosons are interacting or not. It depends on the 
type of the thermodynamic function and the dimensionality of the space. This is in 
agreement with the idea of universality. 
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